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Analysis of the probability distribution of photocount number 
of the onemode stochastic radiation 

Virchenko Yu.P., Vitokhina N.N. 
1. Introduction 

For the description of the registration process of low-intensive electro- 
magnetic radiation by counters of photons, it is necessary to take into 
account quantum effects. Thus, even in the case when the electromagnetic 
field does not contain any noise component, i.e. its amplitude represents 
^ ■ the pure quantum state, the number n of registered photons is random and 

Ph ; there exists the problem of determination of the probability distribution of 

this random variable. This problem becomes very complicated if the regis- 
tered electromagnetic field is the mixture of deterministic field and a noise. 

m ■ 

In this case, its state is statistically mixed and, therefore, it is necessary 
(C~j ; to use the appropriate density matrix. From theoretical point of view, the 

complexity of the description of such a state consists of the adequate choice 
of mathematical model of the electromagnetic noise and of the construction 
on its basis the appropriate density matrix. Simplification of this problem 
arises at sufficiently large typical frequencies of the electromagnetic field. 
In this case, one may use the quasi-classical approximation as it is shown 
in [1]. Therefore, probabilities of registered photon numbers are calculated 
on the basis of the classical (not quantum) probability distribution 

>< : P n = Pr{n = n} = -E>exp[-J], (1) 

1 Til 

which represents the so-called composite Poisson distribution. It is referred 
to the Mandel distribution in quantum optics [2]. Here, J is the random 
variable representing the energy of electromagnetic field absorbed during 
the registration time T. 

Let us consider the model of the quantum photocounter [1] of the 
onemode electromagnetic radiation being completely noise. In this case, 
the appropriate model of the electromagnetic noise is the complex Ornstein- 
Uhlenbeck process as it is proposed in the photodetection theory. Then 
the random variable J is represented by the formula [1] 

T 

J = J[(}= [ ((s) 2 ds 
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where £(s) = £(s) + ifj(s), s G R are trajectories of the complex process 
connected with real Ornstein-Uhlenbeck's processes £ = £ R}, 

fj = {r](t);t G R} being stochastically equivalent and independent. From 
the physical point of view, they correspond accordingly to electric and 
magnetic constituents of the noise electromagnetic field. 

Ornstein-Uhlenbeck's processes are markovian and gaussian and they 
are completely characterized by these properties and their stationary con- 
dition. This class of processes is parametrizated by two numbers v > 0, 
a > 0. Each Ornstein-Uhlenbeck process is completely determined by the 
following formula of the conditional probability density w(xo, to\x, t) of the 
transition from the point Xq G R at to G R to the point x G R at t G R. It 
depends on parameters o and has the following form [3] 



w ^ ') = (i _ e -^-„|) j 6XP {- ~ a (1 - e-2.K-.ol) 

' (2) 

Thus, the one-point distribution density x G R of the process is 

determined by the formula 



/ v \ 1/2 / 

«;(#)= lim w;(xo, tok, ^) = ( — ) exp - 

to^-oo \7T<7/ \ 
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(7 



(3) 



The characteristic function Q(—i\), A G R of the random variable J[£] 
of the process £ is given by the known Ziegert formula [4] , 

Q £ -(A) = Eexp(-AJ[ f ]) = + , )2exp(rT) _ (r _ , )2 exp( _ rT) j . 

(4) 

where r = \Jv 2 + 2A<7. 

Since processes {£(£)}, are independent and equivalent, the gener- 

ating function of the random variable J[Q is found on the basis of equalities 

Q(A) = Q f (A)(^(A) = Qp) . 

From here, we see that the Mandel distribution P n determined by the 
formula (1) and by the probability distribution of the random variable J 
induced by the probability distribution of the process C = {((t)]t G R} are 
very complex. It is easy to obtain the asymptotic formula of the probability 
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distribution P n at T — > 0. It has the form P n /P^ — > 1 where Pi°^ is the 
following Poisson distribution [1] 




exp 



v 



The purpose of the present work is the construction of the convenient 
calculation algorithm giving probabilities P n with the guaranteed accuracy 
at the sufficiently small value of T on the basis of formulas (l)-(4). Our 
approximations are based on the simple idea of the series expansion of the 
Mandel distribution, i.e. we represent it in the form 



oo 



1=0 

It is easy to show that each moment EJ n , n G N is proportional to T n at 
T — > 0. Then one may expect that just such an expansion is appropriate 
for the solution of the above-mentioned problem at small values of T. On 
this way, we should solve two questions. The first, it is necessary to give 
the general foreseeable formula of moments E J n depending on the number 
n as on the parameter. It appears that it is rather problematic to obtain 
such a formula. Instead of this, it is possible, however, to point out the 
specific algorithm of the consecutive calculation of these values. One may 
think that this algorithm is the solution of the first problem. The second 
question consists of the estimation of the series remainder 

oo 1 

£ n Ejn+ ' (6) 

l=N+l 



for any N G N. The presence of such an estimation permits to give the 
guaranteed accuracy of the approximation P^ based on the account only 
of first N terms in the series (5). Below, we solve these two problems. 

2. Sequence algebra of power series coefficients 

Let us introduce into consideration the linear manifold £/ = {f(z)} of 
functions depending on the variable z G C. Each of them is analytic in 
the appropriate circle having nonzero radius and the centre in the point 0. 
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Thus, each function of £/ is presented as the series 

oo 

/(z)=£a»z* (7) 

fc=0 

and each such a series is completely determined by the sequence of coeffi- 
cients (a k G C; k G N+) which is regarded as the infinite ordered collection. 
In this connection, we shall consider also the linear manifold £ = {A} 
of coefficient sequences A = (a k ;k G N + ) corresponding to serieses (7). 
We shall name elements of such sequences as components also. Linear op- 
erations on manifold £ are introduced by the natural way. Namely, let 
A = (a k ; k G N + ) and B = (b k ; k G N + ) be two sequences. Then sequences 
A + B and XA for any A G C are determined by formulas 

A + B = (a k + b k ; k G N+> XA = (Xa k ; k G N+> . 

There exists the natural one-to-one correspondence between introduced 
manifolds £ <-> £/ determined by the formula (7). We designate the 
mapping £ i— > £/ generated by this correspondence using the letter F. 
Further, let us designate the power series of f(z) being the image of the 
element A G £ under the mapping F : £ i— > £/ by F[z|A|. The mapping 
F : £ i— > £^ is obviously linear, i.e. following relations take place 

F[z\A + B} = F[z\A] + F[z\B] , F[z\XA] = XF[z\A] . 

For each sequence A G £, it is possible to consider each of its component 
a k , k G N + as the appropriate projection of the infinite ordered collection 
A. We shall write down this fact by means a k = (A) k , k G N+. 

Let us introduce the binary commutative operation on £ which we shall 
name the convolution of sequences. We designate it by the symbol o. This 
operation is determined as follows. With any pair of sequences A, B G £, 
we associate the sequence A o B having components 

k 

{A oB) k = Y^ djh-j k = 0, 1, 2... . 

3=0 

It is easy to see that the associativity property is fulfilled for the introduced 
convolution operation applied for any three elements A, B,C G £, 

(A o B) o C = A o (B o C) . 
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It concerns also the distributivity property relatively the addition, 

{A + B)oC = AoC + BoC. 

Besides, the following relation 

(XA) oB = X{AoB) 

takes place for any pair A, B E £ and for any A G C. Thus, these equali- 
ties together with the commutative property of the convolution operation 
permits to conclude that the linear manifold £ equipped with the "multi- 
plication" o turns into the commutative algebra which we shall designate 
by the same symbol £. 

In the algebra £, there exists the unity E which is represented by the 
collection E = (1,0,0...) since for any A G £, the following relation takes 
place 

A o E = E o A = A . 

Let us notice that the subalgebra £o = {^4 : (A)q = 0} of the algebra 
£ is its ideal, i.e. it takes place A o B G £o for any A G £o and for any 
B G £. 

There exists the inverse element A^ 1 for any element A £q which has 
the property formulated by the following way 

A- 1 oA = AoA- 1 = E. 
One may find consecutively each its component by the equality system 

(^)o = « 1 

n-l 

(A-% = -clq 1 J2( A ^)n-k a k n = 1, 2, ... . 

k=l 

Let us introduce the following reduced designation of powers of any 
element A G £. We shall write 

A o = E, Al = A, AoA = Al ... , A n o A = A n Q +1 , n = 0,1,2.... 

Each component of the nth power of the element A = (a m ; m G N + ) is 
obtained by the formula 

K) m = fl *i- fl *» m = l,2,.... (8) 

k 1 ...,k n >0 
ki+... + k n =m 
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Thus, it is obvious that if A G £o then the formula 

W) m = E °il- fl in ( 9 ) 

n>jl--iJn>l, 
Jl+...+jn=m 

takes place. Therefore, the relation 

(A n ) m = (10) 

is fulfilled for all m > n. 

The functional F[z|A] is multiplicative relatively the introduced multi- 
plication operation, i.e. the identity 

V[z\AoB] = V[z\A]V[z\B] (11) 

takes place for any pair of collections A and B. It is ascertained by the 
following transformations 

oo 

F[z\AoB] = Y,z m (AoB) m = 

m=0 

oo m oo oo 

= E Z ™ E a ^^fc = E ^ E ^'^m-k = 
m=0 k=Q k=Q m=k 

= (f>^) (f> m M = F[*|A]F[*|5]. 

\fc=0 / \m=0 / 

They are correct in the common convergence area of both series F[z|A] and 
F[A\B]. 

On the algebra £, it is possible to consider "analytical" functions f Q (A). 
Each of them is defined for every A G £ by means of the corresponding 
"power" series 

oo 

fo{A) = Y J C n A n Q c n GC, n G N + . (12) 

Naturally, it has the sense under the condition of the componentwise con- 
vergence. In view of Eq.(10), such serieses are finite for each individual 
component of the function value f Q (A) if the element A is chosen in £q, i.e. 



(/o(^))„ = E c »(^), 



fc=0 



Therefore, the series (12) converges by the definition. In particular, we 
introduce the following function 

oo 

(E-A)? = jT,AS, (13) 

by means of the formula (12), which is defined at A G £o- At last, for any 
function f Q on elements of the algebra £, the following formula is correct 

F[z\f (A)] = f(F[z\A}). (14) 

It is due to linear and multiplicative properties of the functional f[z\-]. 
Here, the function / Q (-) on £ in the lefthand side is determined by means 
of Eq.(12) and the corresponding analytic function on Z determined by 
Eq.(7) in the righthand side is designated by /(•). 

3. Moments of the random variable J 

At first, let us consider the problem about the power expansion on A of 
the generating function Q(X) = Ee~ AJ of the random variable 



J = J CM 



This function is evaluated by the formula 



dt 



n(x) = 4rz/exp(z/T) 

WK ) (r + vf exp(rT) - (r - vf exp(-rT) 1 ; 



on the basis of Eq.(4). Here, r = \/u 2 + 2Xa = vq, 

1/2 



(16) 



We represent the function Q(X) in the form 

Q(A) = e^-^A) 

where 

G(A) = Q-\X)e T = 1 [(1 + qfe" T - (« - ife^] 



Here, we have introduced the "dimensionless" parameter vT instead of T 
which we designated hereafter by the same letter T if it will not cause a 
misunderstanding. 

At first, let us decompose the function G(X) in the series on q powers. 
We have 



G(X) = (iq)- 1 



00 I r ,T\n 00 

d + d'E^r-fo-D'B- 1 ) 



n=0 



n=0 



= (4«) - 1 



Thus, 



d + « 2 )E^r< 1 - (- 1 )") + 2 «E ( J + 



n=0 



= (2«) 



-1 



00 



n=0 

2n+l 



n! 



(«T) 



2n 



n=0 



(2n+ 1)! 



n=0 



(2n)! 



= E 



(9T) 



2n 



n=0 



(2n)! 



1 + 



T 



= E 



T 2 



^(2n) 
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2n 



n=0 



(2n+l) 
T 



1 + 



(2n+l) 



-5 



(17) 



Now, let us substitute the expression (16) of the variable q = (1 + z) 1 / 2 , 
z = 2\(j/v 2 into the formula (17), 



T 2 



«.=() 



1 + (2^TI)( 1 + I 



Then we use the binomial formula. As a result, we obtain 

T 



gw = E 



n=0 



(2n) 



E 

m=0 



n 



1+ (2^T)( 1 + I 



00 z m 00 

= E^™ + « + 2XT 



00 



m=0 



2 * — ' rn 

m=0 



= E 



m=0 



m! 



m 



where 



u 



rn 



= E 



x T 2 » n! 



n=m 



(2n)\ (n — m)\ 



m G N + , 



(18) 



u 



rn 



^ rp2n+l 



nl 



m G N 



+ 



(19) 



^ (2n + 1)! (n - m)\ 

,i=m v 7 x 7 

Consider the sequence W = (w m ; m G N) with components wo = and 

1 ,~T (n . , _ /o , „. A _ T\T (20) 



w m = — e (w m + mv m -i/2 + u TO ) , raGN. 



One may note that e T (wo + ^o) = 1- Therefore, the expansion of the 
function G(X) in the series on A powers is presented in the following form 



G(X) = e 



T 



oo 



z m w, 



m=l 



= e T (l + F[z|iy]) 



2A(j 



2 = 



(21) 



Since (W)o = 0, the sequence W is the element of the algebra £q. 
Therefore, having written down the decomposition (21) in the form 



G(X) = e T [l-f[z\-W]) 



(22) 



and using firstly Eq.(22) and secondly Eq.(13), we obtain the following 
expansion of the function G~ 1 (X) on A powers, 

G-\X) = e- T (1 - F[z\ - W\y l = e- T F [z\X] , 

where X = (E + W)~ l = (x n ; n G N + ), x = 1, 



x 



n 



= yj(-i)' (»-;')„, neN 



1=1 



Thus, we find 



Q(A) = e T G~\X) = t[z\X] 



(23) 



(24) 



substituting the expansion (24) in the expression of the function Q(X). 
Formula (24), in the combination with definition of the variable z, permits 
to us to be convinced that the following statement is valid. 

Theorem 1. Moments M n , n G N of the random variable J are presented 
by the formula 

M n = (-l) n n\(^) x n , neN (25) 



where coefficients Xfi are determined by formulas (18) - (20), (23). 
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□ The formula (25) follows directly from the definition of moments M n 
on the basis of their generated function Q(X). Let us use the expansion 
(24) of Q(X). Then, for any n G N, we have 

*-(-r(?g?>) - 



= (-1)* 



z=0 



We note that each coefficient x n has the sign (— l) n , n = 1, 2, ... due to 
the positivity of moments M n . 

4. Estimations of moments 

For the solution of the accuracy estimation problem of Mandel's distri- 
bution approximations, it is necessary to find some a priori estimations 
of moments M n of the random variable J. This section is devoted to the 
obtaining of such estimations. We divide the obtaining in some simple 
steps. 

Lemma 1. At aG [0, 1/2], the following inequality takes place 

ln(l-a)>-2a. (26) 

□ At a = 0, Eq.(26) turns into the exact equality. At a G [0, 1/2], the 
inequality 

1 

> -2 



1 — a 

of derivatives of Eq.(26) both parts is valid. We obtain Eq.(26) integrating 
last inequality using the equality condition at a = 0. ■ 

L e m m a 2. The inequality 



(2n)! ~ 2 2n {n\y 
takes place. 

□ Following identical transformations are valid 



(2n)! = 2 n n\{2n - 1)!! = 2 2n {n\f f[ fl - ^) = 

i=i ^ 21 ' 
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= 2 2n (n!) 2 exp (J^ ln(l - (2/)" 1 ) j . (28) 

Let us use the inequality (26) for the estimation of the logarithm from 
below, 

ln(l - (2/)- 1 ) > -r 1 , I = 1, 2, n . (29) 



We give upper estimation of the expression l/l considering it as the 
integral sum of the function if (a) = a 



i=\ 
,-i 



n 

Y,) <1 + y ^ = l + lnn. (30) 
z=i i 

The inequality follows from the low estimation of the integral in the right- 
hand side by the rectangular method taking into account that the function 
(f (a) is decreasing at a > 0. 

Now, applying inequalities (29) and (30) for the low estimation of the 
righthand side of Eq.(28), we find 

(2n)!>=-^-. ■ 

en 



L e m m a 3. The following upper estimation of coefficients u m , m = 
1,2,... is valid, 

{T/2) 2m 



(to — 1)1 



where 



oo 



(T/2) 
(n!) 2 

n=0 v ' 



2n 



is the zero order Bessel function on the imaginary variable. 

□ Using the inequality (27), from the definition formula (18), we have 
for any m 6 N that 



- nl - (T/2) 2n , 

^ (2n)! [n - to)! ~ ^ (n - l)!(n - to)! 
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(r/2) 2m ^ (r/ 2 )2(n-m) (r/2) 2m 

- 6 (m-l)! ^ [(n-m)!] 2 ~ 6 (m-l)! °^ 



n=m 



We apply here the elementary inequality (n — 1)! > (to — l)!(n — to)! which 
takes place at n > m. ■ 

Corollary. From Eq.(31) it follows the upper estimation 

Vm < 7r^—ru m = — } (T/2) 2m+1 I (T) (32) 
2m +1 to! 

of coefficients v m , m = 1, 2, ... smce (2/ + 1)! > (2m + 1)(2/)! for all I > m . 

Summarizing, it is possible to assert on the basis of Eqs.(20), (31), (32) 
that the following lemma takes place. 

L e m m a 4. For coefficients w m , the following estimation is valid 

"*» < U e rn ( T / 2 > 2m 'oiTMm) , (33) 
to! (to — 1J! 



(^(to) = 



T TO 

1 + 2^ + r 



L e m m a 5. For modules of coefficients x m , at to > 1, following esti- 
mations take place 

where the function ip(T) is determined by the equality 

ip(T) = Th{2) + y (/ (2) - 1) + 7 (2) (35) 

and 

dT ^ (T/2) 27 "- 1 



^to!(to-1)! 

771=1 V 7 

is the Bessel fubction of first order. 

□ At to > 1, the formula (23) rewrites as follows 

m— 1 /c 

= -w m + ^(-l) fc _1 ^ w m _ Zl _..._ Zfc . (36) 

fc=l h,...,i k >o j=i 

ll + ...+lf^<m 
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Then, using the estimation (33), we have 

m—1 k 

\x m \ < \w m \ + ^ ^2 \ W m-h-...-l k I]IKI< 
fc=l h,...,i k >o j=\ 

l\ + ...+lf.<m 

k=i h,-,i k >o j=i 3 w / 

1 j 1 / 1 

k 

A' 1 \ I 1 

According to the definition of Bessel's functions on the imaginary vari- 
able, the sum in brackets in the last expression is equal to 



oo 



^Z!(^T)! =/l(2) + T(/o(2) " 1)/2 + /o(2)/T= T 

Then, taking in mind of the above obtained estimation and the inequality 
eip(T) > 1, we find 

\x m \ < (r/2) 2 -f>/ (r)) fc (^(T)/T) k = .ff y = (er^(T)/o(T)/4) 

eip{T)I {T) - T 

Corollary. Basing on Theorem 1, from obtained estimations (34), 
we find following estimations of moments M n 



n 



M„ = n!^j \ Xn \< eMT)Io(T) _ T n^ — 2 j . 

R e m a r k. At T — > 0, following asymptotic formulas take place 
u m = [T 2m /(2m)\] (1 + o(l)) , v m = [T 2m+1 /(2m + 1)!] (1 + o(l)) 
Then 

= [T 2m - 1 /2{m - l)!(2m - 1)!] (1 + o(l)) 



w m 
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and, according to the formula (36), the basic contribution in x m gives the 
term with k = 1 at T 0, i.e. x m = [(-l) m (T/2) m ] (1 + o(l)). Therefore, 
since ^(0) = const, the obtained estimation (34) is asymptotically exact 
at T -> 0. 

5. Approximations of the Mandel distribution 
and estimations of their accuracy 

Let us consider the Mandel distribution (1). We present it as the ex- 
pansion on moments 



P - ipy^V- 1 y(-!) 



n! ^ /! n!^ /! 

1=0 1=0 



■M; 



l+n 



or, on the basis of the representation of moments (25), it gives 

(-1)"^ (n + l)\ (2a x 71+1 

1=0 



We determine the sequence of approximations P^ N \ N = 1,2,... of the 
probability distribution which present it with the accuracy up to the iVth 
power of the parameter (aT/v 2 ) (see Remark of Lemma 5) when n < N, 

^^E^(S)" + -. m 

i=o ■ \ / 

Further, we estimate the deviation of the (N — l)th approximation from 
the exact distribution (1). On the basis Eq.(37), we have 



OO 

|p _p(N-l)\ < 1 V 



[n + l)\ (2a 



n+l 



l=N-n x 



1 oo 

= -Y 



l\ 
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N 



Then we use the estimation (34), 

00 71 /r> \ i 



| n n 1 " n! e^(T)-TMl-n)! 2 2 ' ' 1 } 



l=N 
14 



where the variable of the function i/j hereinafter is not pointed out explicitly. 
Let us give separately the estimation of the series 

oo , 

for positive values of £ and for n > N. The following transformations are 
valid 

R ^ ») = c E tt^tC'-" = E C = c c 



<> 2^\iJ(N -OP (i-C) 1+ "-' i-Ct? VJU-C 



TV n /j\r\ / a \ n ^ 

< 



/=0 



Applying the obtained estimation to the inequality (39), we come to the 
following statement. 

Theorem2. The Nth approximation of the Mandel distribution 
presents this distribution with the guaranteed accuracy determined by the 
inequality 

ip _ pw, < ( ;W ^ A2C) W+1 \ 

^ B " 1 - \eij>I (T) -t)\\-2Q ) 

if 

errTj.- 1 
C " ~2^ MT) < 2 ' 



(40) 



In order to obtain the effective algorithm of consecutive calculation of 
Pn N \ it is necessary to solve the last problem. It consists of the creation 
of the method which permits to find consecutively all components of the 
sequence X. 

L e m m a 6. The following formula takes place 

— exp(±a 1 / 2 T)J = m\e ±T Ri(T) (41) 
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where polynoms R^(T) of the m degree on the variable T are determined 
by the recurrent relation 

R - +i{T) = ± n^)f^ { - 1)l My R ^-' (T) ' meN+ (42) 

and by the condition at m = 

Ro(T) = l. (43) 

□ The proof is realized by the induction on m. At m = 0, we obtain 
(42). We shall build the induction step. Let Eq.(41) takes place at the given 
value m G N + . Then we shall introduce the function R m+ i(T) according 
to 

exp (ia^T 



a=l 



= ± 



da 7 



1 exp(±a 1 / 2 T) 



2a 1 /2 



a=l 



Using the induction assumption and the formula 



(«" 1/2 ) (i) = (-i)'||«- ,2i+1,/2 . 

we accomplish the ra-fold differentiation, 
(m + l)\e ±T R m+1 (T) = 

=4g(T)(£-(— ))„(£ 



a=l 



- ±mn — 2^ 7j R m -i\i )^qj 

1=0 



Whence formulas (41), (42) follow and although the fact that the function 
R m+ i(T) is the polynom of the m + 1 degree is proved. ■ 

Theorem 3. Components of sequences U = (u m ; m G N) and 
y = (v m ', m G N) are represented in the form 

um = \m\ (e T R+(T) + e~ T R m {T)) , (44) 
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v m = i(m + 1)! {e T R + m+l {T) + e~ T R- +1 (T)) 

where polynoms i? m (T) are calculated by the formula (42). 

□ Using Eq.(41), from the definition of coefficients u m , we have 



(45) 



u 



rn 



n=m 



rp2 



n 



Qm a n T 2n 

~doF l ^ (2n)\ 



d 



rn 



da 



in 



ch{a l ' 2 T) 



(2n)! (n — m)\ 



n=m 



a=l 



a=l 
1 



2\da 



in 



exp(a 1/2 T) + exp(-a x 




a=l 



On the basis of the formula (41), we have although 



da 



rn 



TO 



m 
J, 

TO 
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m—l 



da 



rn- 



exp (±a 1/2 r) 



a=l 



a=l \ dal 



a 



-1/2 



a=l 



= mI e^£(-i)<jj£_ ( (D = ±|(m + l)!e^ +1 (T) 



;=o 



2 2 *(Z!) i 



T 



Further, on the basis of the obtained formula, from the definition of coef- 
ficients v m , we have 



oo 



V 



in 



T 



2n+l 



ni 



(2n + 1)! (n — m)! \ <9ai 



5 to ^ a «T 2n+1 



(2n + 1)! 



a=l 



d m sh(a 1 / 2 T)' 



1 / d r 



da 7 



1/2 



a=l 



2 Vda r 



exp (a 

1/2 T ) 

exp(— a 

i/2 T y 



a 



1/2 



a 



1/2 



J / a=l 



1 



= -(m + 1)! (e T R+ +1 (T) + e- T i? m+1 (T)) 



Corollary. Taking into account Eqs.(44),(45) and Eq.(20), compo- 
nents of the sequence W are represented as follows 

»» = \R + JT) (i + ^) + + 
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+ e- 2T 



1^ /rriS ( m\ m + 1 



m G N . (46) 



6. Approximate formulas 

Thus, Theorem 3 permits to calculate successively all approximations 
P^ of probabilities P n using formulas (46). In this section, we shall give 
some first approximations of the probability distribution P n by means of 
this algorithm. 

At first, we find the explicit form of polynoms R^(T), m = 0, 1, 2, 3, 4, 

flf(T) = l, Rf(T) = ±T/2, 

R±(T) = ±(T/2 2 ) (R±(T) - 1/2) = (T/2 3 )(T T 1), 
Rf{T) = ±(T/2-3) (i^(T) - i?f (T)/2 + 3/8) = ±(T/2 4 -3) (T 2 t 3T + 3) , 
flf(T) = ±(T/2 3 ) (i?3 (T) - i£(T)/2 + 3i?f (T)/8 - 5/16) = 
= (T/2 7 • 3) (T 3 t 6T 2 + 15T t 15) . 
On their basis, we calculate components it> m , m = 1, 2, 3. 

^i = f, W2 = ^(2r 2 -2r + i- e - 2r ) , (47) 

w 3 = _L (2T 3 - 6T 2 + 9T - 6 + 3(T + 2)e~ 2T ) . (48) 
At last, we find components x m , m — 1,2,3, by the formula (36), 

#0=1? X\ = —W\ , X2 = — W2 + w\ , £3 = — 1^3 + 2w\W2 ~ U)\ 

or after the substitution (46), (47), we have 

^i = -|, x 2 = i (2T 2 + 2T - 1 + e- 2T ) , 

x 3 = -—^5 [2T 3 + 6T 2 + 3T - 6 + 3(3T + 2)e~ 2T ] . 

Explicit expressions of components x m , m = 0,1,2,3 permit to us to 
write the approximated expressions of n-photon registration probabilities 
for n = 0, 1, 2, 3 up to the third order (N = 0, 1, 2, 3). All probabilities P^ 
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of some higher values n > 3 are equal to zero in this approximation. We 
result consequtively those expressions for each order of the approximation. 

For N = 0, we have that only the probability Pq ^ = 1 differs from zero. 

At N = 1, we obtain 

P (i)_ v f—\ Xl -i-— p^-ZL 



Z=0,1 



At A = 2, we have correspondingly 

^ = £ (I)'— l-S+(^) 2 ^ + 2T-l + e -), 

i=0,l,2 V J 



P? = ^-2{^)\2T> + 2T-l + e^), 

pf = (-^) 2 (2r 2 + 2T-l + e- 2r ). 
At last, for N = 3, we get following approximate formulas 



f> = £ 

Z=0,l,2,3 



Xi = 



= 1 - 



oT ( O \ 2 



IS 2 



+ (^) (2T 2 + 2T - 1 + e- 2T ) - 

- - (^2 J [2T 3 + 6T 2 + 3T - 6 + 3(3T + 2)e" 2T ] , 

^^-Ec+Dfl)'"— = 

Z=0.1.2 ^ ' 



=3-- 2 (2^) 2 ( 2T2+2T - i+e " T ) + 

+ 2 (-^) [2T 3 + 6T 2 + 3T - 6 + 3(3T + 2)e~ 2T ] , 

Z=0,1 ^ ' 



Z=0,1 

2 



= (^)V 2 + 2T-l + e--)- 



19 



- 2 [2T 3 + 6T 2 + 3T - 6 + 3(3T + 2)e" 2T ] , 

P 3 (3) = ^ (^2) 3 + 6T 2 + 3T - 6 + 3(3T + 2)e~ 2T ] . 

7. Conclusion 

It is necessary to point out some lacks of the problem setting and its 
solution in the given work. 

Our research is based on the expansion (5) which is directly connected 
with the expansion of the generating function (15) on A powers. The 
convergence radius of last expansion is equal to the convergence radius 
r(T, z/, a) of the power series on A of the function Q(X) which is meromor- 
phic only. It has poles in the complex plane A G C. It means that the 
distance to nearest of these poles is equal to r(T, z/, a). It is clear that 
r(T, z/, a) — > when parameters T or a increasing. At the same time, such 
restriction of the convergence domain of the expansion (5) is not caused by 
a physical reason, i.e. the divergence of the series (5) is not connected with 
the presence of anything qualitative change of the probability distribution 
P n behavior when parameters T, z/, a varying. Therefore, it is desirable to 
get rid of this lack of the expansion (5). 

Other lack is connected with the above mentioned one. With irreversibil- 
ity, we make a mistake when calculating the convergence radius by the esti- 
mation of the residual series (6) basing on a priori estimations of moments 
EJ TO , m G N. Actually, only its low bound is obtained by this way and it 
is sufficiently coarse as one may see from our work. 

It is possible to consider mentioned lacks as ones generated by our 
method when the problem solving. The following lack concern with the 
unsuccessful problem setting. 

Let us pay attention that obtained formulas permit to calculate the 
registration probability of everyone concrete photon number n. However, 
they do not permit to calculate (to estimate) the registration probability 
in the case when this number is indefinite, i.e. in the case when it is 
not known exactly and, therefore, the number n is the parameter of the 
problem. It is connected with the fact that their analytical representation 
becomes very tedious when the approximation order increasing. 
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